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Magnetic field applied to the quantum dot coupled between one metallic and one superconducting 
electrode can produce a similar effect as has been experimentally observed by Meservey, Tedrow and 
Fulde [Phys. Rev. Lett. 25, 1270 (1970)] for the planar normal metal - superconductor junctions. 
We investigate the tunneling current and show that indeed the square root singularities of differential 
conductance exhibit the Zeeman splitting near the gap edge features V = ±A/e. Since magnetic 
field affects also the in-gap states of quantum dot it furthermore imposes a hyperfine structure on 
the anomalous (subgap) Andreev current which has a crucial importance for a signature of the 
Kondo resonance. 



I. INTRODUCTION 

Already in early days of the tunneling spectroscopy 
it has been shown that magnetic field B (which couples 
to spin of the charge carriers) is in superconductors re- 
sponsible for splitting the square root singularities of the 
tunneling conductance pj] by the Zeeman energy 2^bB, 
where /xs is the Bohr magneton. This Meservey-Tedrow- 
Fulde (MTF) effect has been observed experimentally in 
the thin superconducting aluminum films applying par- 
allel magnetic field so that orbital diamagnetic effects 
could be avoided. Similar qualitative results have been 
recently noticed in the measurements of c-axis tunneling 
for the layered high temperature superconducting com- 
pounds 

We argue that the MTF effect should be also feasi- 
ble in various nanostructures consisting of a quantum 
dot (QD) placed between one metallic and one supercon- 
ducting electrode. Zero-dimensional character of QDs in 
a natural way eliminates the influence of orbital effects 
therefore magnetic field would affect the charge transport 
only through the Zeeman term. This can in turn manifest 
itself in the differential conductance. Roughly speaking, 
the charge current flows if an external bias V exceeds the 
energy gap A (necessary to break the Cooper pairs into 
individual electrons) thereof the resulting conductance 
has a low voltage onset near the gap edges eV = ±A. In 
presence of a magnetic field these gap edge singularities 
are going to split (see section III). 

More detailed analysis of the charge tunneling [3] in- 
volves however also the additional (anomalous) channels 
due to mixing of the particle and hole excitations in su- 
perconductors. In particular, even at subgap voltages 
\eV\ < A the mechanism of Andreev reflections provides 
a finite contribution to the conductance. Since the An- 
dreev mechanism is very sensitive to location of the in- 
gapQD statesjj, Ja [SiHJa ^.nd the on-dot correlations 
Em HH, M, M, MmKM, [13, m we shaU explore the 
influence of magnetic field on such subgap conductance. 
In section IV we discuss a hyperfine structure for the 
Andreev conductance neglecting the correlations. In the 



next section V we extend our study taking into account 
a finite value of the on-dot repulsion U. We show that 
appearance of the low temperature Kondo resonance en- 
hances the zero bias conductance and this feature under- 
goes the Zeeman splitting when magnetic field is applied. 

As concerns some practical aspects, there have been 
considered the proposals for using the magnetic field 
tuned Andreev scattering as an efficient cooling mech- 
anism in two dimensional electron gas - superconductor 
nanostructure (Toj . There is also considered a possibility 
to use the, so called, Andreev quantum dot as a magnetic 
flux detector [20j . 

II. THE MODEL 

For a general description of transport phenomena 
through a nanoscopic island placed between external 
leads one should consider a quantized multilevel struc- 
ture of QD [Hj. However, in the case when a level spac- 
ing is smaller in comparison to QD hybridization with 
the electrodes one can restrict to a simplifled picture of 
the Anderson model [i[Ei,[il[ni 

CJ 

+ E E dUua,3 + V^p cL_^d,) . (1) 

k,(T I3=N,S 

Operators (dj.) denote the annihilation (creation) of 
electron whose energy level is ed,a and U is the on- 
dot Coulomb repulsion between opposite spin electrons. 
The last terms describe hybridization of QD with the 
normal (/3 = N) and superconducting {j3 = S) elec- 
trodes. Magnetic fleld eventually shifts the QD level by 
£d,a = £d — gcr^J'BB, whcrc the spin-dependent coefficients 
are defined as — 1 and = — 1. 

Hamiltonian of the normal (metallic) lead is taken 

as Hn = Sk.cr CkJvCkcrArCko-JV whereas for the super- 
conducting electrode we choose the usual BCS form 

^s = Ek,^^£scL5Cko-s-Ek (^^TsC-kiS + with 
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an isotropic energy gap A. The relative energies = 
(sk/3— 3<tA*-B-S)~M/3 s-re measured from the chemical poten- 
tials fi/s. We shall focus on the wide band limit |Vk/3| ^I? 
(where —D < eu/s < D) and consider a small exter- 
nal voltage V, which detunes the chemical potentials by 
fJ-N^fJ-s = eV inducing the charge flow through N-QD- 
S junction. We assume \eV\ to be much smaller than 
level spacings typical for the realistic QDs [21J so that 
applicability of the model H]) can be justified. 

Let us start by establishing the QD Green's function 
in the equilibrium situation, i.e. for V = 0. Fourier trans- 
form of the retarded Green's functions can be formally 
expressed by the Dyson equation 
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where S° ^ denotes the selfenergy of noninteracting QD 

{U — 0) and S^,^ accounts for the correlation effects. For 
a simple understanding of the MTF effect it would be 
helpful to focus first on the uncorrelated QD when the 
selfenergy is known exactly. Further corrections due to 
S^g. contribute a renormalization of the spectral func- 
tion [3] whose impact on the charge transport will be 
discussed separately in section V. 

For convenience we introduce the hybridization cou- 
pling Tp = 27r^j^ |Vk/3p5(w — Ek/s) and define the fol- 
lowing spin-dependent energy oj^ = to + ga-nsB. Imagi- 
nary part of the selfenergy ^ for po-l < A is given by 
ImS° ^.(cj) = -i^l while at large energies D>\u}c\>A 
it takes the following form [l^, [l^ 
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The corresponding real parts can be determined using 
the Kramers-Kronig relations. 

Imaginary part of the selfenergy Yp^ ^ has thus the 
square root singularities at energies w = ±A ± ^irB, so 
in presence of magnetic field there are altogether 4 such 
points. They show up as kinks in the spectral function 
= Y.a Pd,<y{'^)i where 



Pd 



,.H = --im((4;4» 



(4) 



We shall see below that appearance of such characteristic 
points leads to the MTF effect observed in the tunneling 
conductance. 



FIG. 1: The differential conductance G{V) versus bias voltage 
V for N-QD-S junction. Notice a splitting of the gap-edge 
singularities around eV = ±A induced by magnetic field B. 
We used for computations £^ = 0, U — 0, rjv = A, rs = 0.1A, 
T = 0.01 A assuming A = 0.1D. 



here on a situation with the spin sensitive transport due 
to magnetic field. The steady charge current is defined as 



-dt (X]k,(T ^kcrAf^ktr 
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carry out the time derivative and determine the expec- 
tation value using the nonequilibrium Keldysh Green's 
functions. 

In analogy to the standard Blonder-Tinkham-Klawijk 
theory [H we express the current as composed of two 
contributions 



IiV) = h{V) + lA{V). 



(5) 



The first part Ii{V) stands for a contribution which 
at low temperatures appears practically outside the en- 
ergy gap \eV\ > A. Its magnitude is expressed by the 
Landauer-type formula 



dujTi,^{uj) [f{cu + eV)-f{uj)], (6) 



where /(w) = [I + exp{uj/kBT)]~ . The transmittance 
ri,o-('j-') is nonvanishing only outside the energy gap 
jcjcrl > A and is given by the following parts of the re- 
tarded Green's functions H 
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The second part in (0 
of Andreev reflections I 
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I originates from the mechanism 
lilll 



dcuTAA^) [f(oJ + eV)-f{uj-eV)] . (8) 



III. MESERVEY-TEDROW-FULDE EFFECT 



To compute the tunneling current we adopt the formal- 
ism outlined in the previous studies extending it 



Its transmittance is finite even inside the energy gap 



(9) 



3 




FIG. 2: Zeeman splitting of the bound Andreev states for 
the QD located in the center of superconducting gap Sd = 0. 
Upper panel illustrates the density of states pdii^) and the 
bottom figure shows differential conductance of the in-gap 
current. For computations we used FAr^O.lFs, fisB — Q.lTs 
assuming Fg = O.OID and U — 0. 
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FIG. 3: Differential conductance Ga{V) of the in-gap An- 
dreev current as a function of the bias voltage V and the QD 
level Ed- We used for computations Fs = O.OID, Fjv — O.lFs, 
T = O.OlFs and set the magnetic field ^fisB = O.lFs. The 
conductance is expressed in units of Ae^ /h. 



et al [iJI- In such case /i can be completely discarded 
from our analysis. Usin g Q we obtain the selfenergy 
^ simplified to H, [H, [ST 



Physically such process occurs when an incident electron 
from A'^ electrode (of arbitrary energy) is converted into 
a pair on QD (with a simultaneous reflection of a hole) 
and it propagates in S electrode as a Cooper pair. This 
anomalous Andreev current is closely related to the off- 
diagonal order parameter induced in the QD (proximity 
effect) [l3,li3- 

Figure ([T]) illustrates the influence of magnetic field 
on the total differential conductance G{V) — -^I[V) 
obtained for N-QD-S junction. We clearly notice the 
Zeeman splitting of the square root singularities resem- 
bling the former experimental observation for N-I-S (I- 
insulator) junction \V\. However, in a present case the 
conductance does not saturate to a finite value far out- 
side the gap \eV\ ^ A because the QD spectrum spreads 
only nearby Ed (usually in realistic multilevel QDs there 
would be seen the quantum oscillations of G{V) [2l|). 
The in-gap features related to the Andreev current are 
discussed in the next section. 



IV. MAGNETIC FIELD EFFECT ON THE 
ANDREEV CURRENT 

The mechanism of Andreev reflections transmits the 
charge current even for the subgap voltages. To focus 
solely on this anomalous current it is convenient to con- 
sider the extreme limit A — > oo as proposed by Tanaka 



(10) 



Upon neglecting the Coulomb correlations one can an- 
alytically determine the Green's function ([2]), where the 
spin dependent spectral function ([4]) acquires the BCS 
structure [lj| 
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with a quasiparticle energy Ed = ^/e^ -I- (Ts/'^Y- The 
in-gap QD states (often referred as Andreev hound states) 
form around ±£'d±/iB_B as illustrated in the upper panel 
of figured! Their line broadening is given by rAr/2 and 
in absence of magnetic field the particle-hole splitting is 
controlled by Ts [H, [3 (the dashed fine in figure [2]). 
Magnetic field further enforces the Zeeman splitting of 
these in-gap states. 

Above mentioned behavior has an indirect effect on 
the off-diagonal parts of the Green's function ([2]) which 
in turn determine the Andreev transmittance. In the 
limit A ^ CX3 © reduces to 
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FIG. 4; Differential conductance Ga{V) of the in-gap An- 
dreev current as a function of the bias voltage V and magnetic 
field B for the QD level Sd^O and riv = 0.irs, T = O.Oirs . 
Dark areas denote the regions where Ga approaches the value 
4e'^/h. 



The subgap Andreev conductance Ga{V) — -^Ia{V) is 
thus characterized by a four peak structure as shown in 
the bottom panel of figure [2] Obviously the weights of 
particle and hole peaks of the spectral function (jlip as 
well as their weights in the Andreev transmittance (fT2|) 
depend on the QD level Ed- Variation of the Andreev 
conductance with respect to {V, Sd) is plotted in figure 
131 We can notice that optimal conditions for the sub- 
gap current occur when the QD level is located near the 
energy gap center, otherwise the proximity effect is less 
efficient. 

On top of the particle-hole structure seen in the An- 
dreev states there is an additional Zeeman splitting 
brought by magnetic field. In figure |4] we sketch the An- 
dreev conductance in {V ,B) plane for £d = 0, where the 
dark areas correspond to a maximal value Ae^/h. There 
appears a characteristic diamond shape marking the posi- 
tions of such maximal conductance Ga{V, B). We believe 
that this hyperfine structure could be probed experimen- 
tally. 

To complete the discussion of the subgap Andreev cur- 
rent we briefly comment on a possible influence of an 
asymmetry between the hybridization couplings Fjv, Fg. 
We explore for this purpose the zero bias conductance 
GAiy = 0). At low temperature we find from equation 
(ini) that 



Ga{Q) = 
4e2 



(13) 
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In figure [5] we show the influence of magnetic field on 
the zero bias Andreev conductance for several values of 
the asymmetry rate Tn/Ts- If Tn/Ts ^ 1 then a line- 
broadening of the Andreev states diminishes so in con- 
sequence the particle and hole peaks become well sep- 
arated. Under such conditions the subgap conductance 
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FIG. 5: The zero bias differential conductance Ga(0) as a 
function the magnetic field S for several relative values of 
Fiv/Fs. We used for computations ed = assuming T-^O. 



has maxima around the quasiparticle states at ±r5/2 
(where the ideal conductance 4e^//i is reached). Let 
us recall, that in absence of magnetic field the equa- 
tion (fT51) reproduces for = the well known result 

'2r„r^ 



4e" 
h 



Ga(0) = 
Fs^F^ it yields G^(0) 



/ 2ijvr^V 



|. For the symmetric coupling 



V. INFLUENCE OF THE COULOMB 
CORRELATIONS 

In the limit IS. ^ co the selfenergy ^ becomes a 
static quantity (fTU)) therefore the role of superconducting 
lead can be exactly replaced by the on-dot gap parameter 
Ad = Fs/2. Instead of H]) we can thus use the following 
auxiliary Hamiltonian 



U fid^ndi, (14) 



k,(T 

\du\dl 



Addld] + h.c. 



which turns out to be very convenient for investigating 
the correlations. Tanaka and coworkers [il,[l3 were able 
to rigorously prove that the selfenergy must have a 
diagonal structure due to invariance of Uhd^fidi term on 
the Bogoliubov-Valatin transformation. 

In the remaining part of this section we shall focus 
on the subgap Andreev current transmitted through the 
correlated QD. The matrix Green's function ^ simplifies 
in the limit A^oo to the following (exact) structure 



(15) 



TV,- 



Influence of the correlations have been so far analyzed for 
the Hamiltonian ([1]) using various techniques [§, [13, [H!, 
Il2l . [3, [l^ . Here we estimate the diagonal selfenergy 
Sat^ct^) within (|14p by the equation of motion method 
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where ^'^ (lo) are given by [22 
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FIG. 6: Spectral function of the correlated QD obtained for 
Ed ~ — l.SFs, U = lOFs, Fiv = Fs and temperature T = 
10~^Fs (<< Tk) in the limit A ^00. Solid line corresponds 
to fiBB^rs/3. 
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FIG. 7: Effect of magnetic field on the subgap Andreev con- 
ductance Ga{V) in the Kondo regime with Ed = — I.SFat, 
(7 = lOFiv, Fs = 5F]v and T = IQ-^^Fjv < Tk- Notice ap- 
pearance of: (i) the zero-bias Kondo anomaly (showing the 
Zeeman splitting for /^sB = F]v/3, (ii) the quasiparticle peaks 
at |el/|~i5£j, and (iii) Coulomb satellite peaks near \eV\ ~ U. 



Approximation (|16ll7p qualitatively reproduces the 
following properties caused by on-dot correlations: (i) 
the charging effect and (ii) a possible appearance of the 
Kondo resonance for temperatures smaller than Tk = 
^^^^exY>{ned {ed + U) /UTj^}. The latter one is related 
to screening of the quantum dot spin by itinerant elec- 
trons of the metallic lead. In the case when energy level 
Ed is located slightly below /i^r the hybridization VkAr 
induces effectively antiferromagnetic interaction between 
the QD and metallic lead. In consequence the bound 
singlet state can be formed giving rise to the resonance 
at Lo = fiN for temperatures T < Tk- Magnetic field 
eventually splits this resonance as illustrated in figure [51 

Any features present in the QD spectrum are further 
showing up in the measurable differential conductance. 
This is also valid for the Kondo resonance. Since it forms 
near the chemical potential fiN therefore its signatures 
appear predominantly in the low voltage current. In fact, 
it has been shown that Kondo resonance enhances at low 
temperatures the zero bias Andreev conductance [1, [l^ , 
however its magnitude remains much smaller than the 
unitary limit value 2e^//i typical for N-QD-N systems in 
the Kondo regime. In the present context we emphasize 



that magnetic field enforces the Zeeman splitting of the 
zero bias Andreev anomaly in much the same way as it 
affects the zero bias anomaly for the QD coupled to both 
metallic leads [H, [l^ . 

The zero bias enhancement of the Andreev conduc- 
tance is a feature whose presence might be difficult to 



notice [q 
fulfilled 




I15j unless some stringent requirements are 
25j . It turns out that optimal conditions for the 
low temperature enhancement of Ga {V ~ 0) take place 
when Ts is comparable to F^r (see figure [8|) and Ed is 
located slightly below the energy gap center. For an in- 
creasing asymmetry between the hybridizations Fat, F5 
the magnitude of low voltage Andreev conductance di- 
minishes (similarly as we have been shown in section IV 
upon neglecting the correlations). On the other hand, 
for Ed moving far aside from the superconductor's gap 
center the proximity effect becomes weakened and the 
overall Andreev conductance is again suppressed. 

In general it seems that an interplay between the on- 
dot pairing (absorbed from the superconducting elec- 
trode) and the Kondo state (due to screening of QD spin 
by the metallic lead electrons) has the same character as 
a competition of superconductivity versus magnetism in 
the solid state physics. Since this is outside the main 
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Ga(V) 




FIG. 8: The difFerential Andreev conductance Ga{V) (in 
units 4e^//i) as a function of the bias voltage V and the asym- 
metry ratio Ts/Tn- We used the same set of parameters as 
in figure El 



scop e of the present topic we shall discuss it separately 
[25| . A combination of the Kondo physics, superconduc- 
tivity and the Zeeman polarization is a complex problem 
and to our knowledge only few papers have so far at- 
tempted to address this challenging issue 26]. 



VI. SUMMARY 

We have explored the effect of magnetic field on charge 
transport through the quantum dot attached to one nor- 
mal and one superconducting electrode. For a bias volt- 
age F~±A/e we find the Zeeman splitting of the square 
root singularities in the differential conductance. This 
resembles the experimental result of Meservey, Tedrow 
and Fulde observed in the N-I-S junction [l| which for 
the N-QD-S structures it seems rather easy to achieve. 

We have extended our study also on the in-gap An- 
dreev current. Due to the proximity effect the particles 
and holes of the quantum dot get mixed and effectively 



the spectrum acquires the BCS-like structure pT|) . Dif- 
ferential conductance GAiV) of the in-gap current indi- 
rectly probes such structure of the bound Andreev states. 
We have shown that magnetic field leads to appearance 
of four peaks via the combined particle-hole and Zeeman 
splittings. We hope that this result might stimulate a 
search for the experimental detection of above mentioned 
structures. 

Moreover, we have explored influence of the on-dot 
Coulomb interactions on the subgap Andreev current as- 
suming the extreme limit A— >oo. In general, the on-dot 
correlations contribute to the QD spectrum: (i) appear- 
ance of the Coulomb satellite near lj = ed,^ + Sd,i + U 
(charging effect), and (it) at sufficiently low tempera- 
tures can produce the narrow Kondo resonance at the 
chemical potential fi^. Magnetic field imposes the hy- 
perfine splitting onto such spectrum in a similar way as 
has been observed in N-QD-N junctions |2j|. The Kondo 
effect alone is exemplified in the zero bias Andreev con- 
ductance where under appropriate conditions [1^ a low 
temperature enhancement can be seen if F^v Fg and 
the gate voltage tunes nearly to the energy gap center. 

It would be of interest to use some more sophisticated 
methods for treating the on-dot interaction U in order 
to check whether there exist a minimal magnetic field 
necessary for splitting the Kondo peak (as theoretically 
predicted for N-QD-N junctions ^23i|) observable in the 
Andreev conductance. One can also study QD coupled 
with d-wave superconductor, where the square root sin- 
gularities are replaced by weaker kinks. We think that 
the Meservey- Tedrow-Fulde effect would be observable 
there too (but in a less pronounced manner) whereas the 
subgap conductance might qualitatively change. 
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